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Abstract 

Massless scalar electrodynamics is studied at high temperature and zero chemical po- 
tential. I derive the free energy to order A^, Ae^ and e'^ by effective field theory methods. 
The first step consists of the construction of an effective three-dimensional field theory that 
is valid on distance scales 1 /T and whose parameters can be written as power series in 
A and e^. These coupling constants encode the contribution to physical quantities from the 
scale T. The second step is the use of the effective Lagrangian and perturbative calculations 
yield the contributions to the free energy from the scale eT ^ ^/XT . 

PACS number(s): ll.lO.Wx, 12.20.DS 



1 Introduction 



Quantum field theory at finite temperature have received considerably attention in recent 
years for a number of reasons. Two major applications are phase transitions in the early 
Universe and QCD at high temperature. 

It was first observed by Kirzhnitz and Linde that symmetries which are spontaneously 
broken at zero temperature are normally restored at high temperature The phase tran- 
sition which has been studied most extensively is the electroweak phase transition, which 
takes place when the Universe has a temperature of approximately 200 GeV, and where 
the Higgs field acquires a nonzero expectation value. One reason for investigating the 
electroweak phase transition is that it may be responsible for the present baryon asym- 
metry. In order to explain the excess of baryons over antibaryons, the electroweak phase 
transition must be sufficiently strongly first order 0. For realistic values of the Higgs 
mass, this is not the case in the minimal standard model [0, and so one is led to inves- 
tigate extensions, such as two Higgs doublet models |0] and the minimal supersymmetric 
standard model 0]. 

The hadrons in QCD also undergoes a phase transition when the temperature is raised 
to approximately 200 MeV. The hadrons turn into a plasma of free quarks and gluons. 
This quark-gluon plasma may be produced in heavy-ion collisions in the next generation 
of colliders. Hence, the study of QCD in the deconfined phase is very important. 

It is now a well-known fact that bare or naive perturbation theory at high tempera- 
ture breaks down for soft external momenta /c (/c ~ gT, where g is some generic coupling 
constant) due to infrared divergences 0. Leading order results for physical quantities 
such as screening masses and damping rates receive contributions from all orders of per- 
turbation theory. In order to do consistent perturbative calculations, it is necessary to 
apply the resummation program of Braaten and Pisarski or extentions thereof. In this 
approach one uses effective propagators, and in the case of non-Abelian gauge theories 
one uses effective vertices as well. This effective expansion resums the infinite subset of 
diagrams which contributed to physical quantities al leading order, and hence it is truly 
perturbative. 

Resummation methods have been used by several authors to study the symmetry 
restoration at finite temperature of spontaneously broken gauge theories. Hebecker has 
studied the Abelian Higgs model to two- loop order in resummed perturbation theory 
while Fodor and Hebecker have analyzed the electroweak phase transition, also in the 
two-loop approximation [0 . The Abelian Higgs model has been used a testing ground for 
different resummation methods in hot gauge theories [Q. However, it is also of interest 
in its own right, since it provides a model for a relativistic superconductor. The phase 
transition is then a transition from a superconducting phase (Higgs phase) to a normal 
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conducting phase (symmetric phase). 



In some cases, as in the study of static quantities such as the effective potential, Arnold 
and Espinosa have invented a simplified resummation scheme 0. For static quantities, 
it is sufficient to dress the bosonic zero-frequency modes by thermal masses. For all 
other bosonic modes as well the fermionic modes, the Matsubara freqencies provide the 
necessary infrared cutoff. Arnold and Espinosa have used this resummation approach 
to study phase transitions in various gauge theories. It has also been used to compute the 
free energy in gf^^'^-theory to order by Parwani and Singh (the result to order g'^ was 
first obtained by Frenkel, Saa and Taylor |]TT[). Corresponding calculations in QCD have 



been carried out by Arnold and Zhai to order g^ |jT2|, and extended to g^ by Kastening and 



Zhai |T^. Finally, Blaizot et al. have used the method to compute the electric screening 



mass in scalar and spinor electrodynamics to order and e^, respectively [|T4 



In recent years effective field theory methods have invaded the area of high tempera- 
ture field theory. This was initiated by the papers of Ginsparg and Landsman [15,16]. The 
basic idea is as follows. In the imaginary time formalism there is a summation of Matsub- 
ara frequencies in loop diagrams. These frequencies can be viewed as masses proportional 
equal to 27rnT for bosons and {2n + 1)ttT for fermions, where n is positive integer. Hence, 
there is a mass hierarchy, and the fermionic modes as well as the nonzero bosonic modes 
ought to decouple at distance scales R ^ 1/T, analogous to the decouphng of heavy 
particles at zero temperature and low energy according to the Appelquist-Carrazone the- 
flTTl. 



orem 



There are several ways of constructing effective field theories at finite temperature. At 
the one-loop level, one would integrate out the heavy modes. In this context this means 
that the coupling constants are defined by computing the effects of the heavy modes 
on static correlators. Hence, only heavy modes circulate in loop diagrams. Beyond one 
loop, one must be more careful. It has been demonstrated by Jakovac that the above 
procedure generates non-local terms which cannot be expanded in local operators, and so 



it is difficult to obtain £gff [Q. The reason is that there can be light and heavy modes 
simultaneously in multi-loop diagrams, and that one must take this into account in order 
to obtain the effective Lagrangian. 



Kajantie et al. have put forward an effective field theory approach, which resolves 
the above problems. The coupling constants are now determined by a matching procedure. 
One calculates static correlators in the full and in the effective theory, and demand they 
be the same. At the one-loop level, this method coincides with the original approach to 



dimensional reduction |]I6[. Beyond the one-loop approximation, one omits the contribu- 
tions to a correlator where all the modes circulating are light, but contributions involving 
both heavy and ligth particles are retained. The effects from a loop diagram involving 
only static (light) modes will be taken care of by the effective Lagrangian. In many cases. 
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one needs the correlators at zero external momentum. One can then directly apply the 
effective potential, since this is the generator of one-particle irreducible Green's functions 
at zero external momentum. One advantage of using the effective potential is that this 
often is easier to calculate than perorfming a direct calculation of loop corrections to the 
correlators. 



There is an alternative effective field theory approach due to Braaten and Nieto [pO[ 



Instead of explicitly constructing an effective theory for the zero mode one proceeds as 
follows. One writes down the most general Lagrangian in three dimensions, which respects 
the symmetries of the theory and includes all relevant degrees of freedom (i.e. no fermionic 
fields are present in Cqq). The coupling constants are again determined by computing 
static correlators in the two theories and require that they match. One advantage of this 
approach is that we dot not explictly divide loop corrections to a static correlator into 
contributions from light and heavy modes. This method has been applied by Braaten 
and Nieto both to (7^$^-theory and QCD They have confirmed the results of 



calculations of the free energy to order in these theories, first obtained by resummation 
methods by Parwani and Singh |T^, and by Zhai and Kastening [|T3], respectively. They 



also computed the screening mass squared in (7^$^-theory to order g'^, which is a new 
result. The present author has applied the methods to scalar ||2^ and spinor QED 



and reproduced the results by Blaizot et al. ||T^ for the electric screening mass in the 
two theories to order and e^, respectively. Similarly, the free energy result in QED to 
order found by Zhai and Kastening agrees with that obtained by effective field theory 



methods |2^. This effective field theory aproach has also been used to calculate the 
screeening mass squared in g^^'^ theory to order g^ |24 



In contrast with other theories, the literature on SQED at very high temperature 
(T ^ Tc) is somewhat sparse. There are no papers specifically devoted to the calculation 
of the free energy at high temperature, although the result to order and A'^''^ may be 
easily extracted from existing literature ^ . In this work, I apply the effective field theory 
methods of Braaten and Nieto to compute the free energy in SQED to order Ae^, and 
e^. The result then represents the present calculational frontier. 

The plan of the present paper is as follows. In section II, we discuss dimensional 
reduction and how effective field theory methods can be applied high temperature ther- 
modynamics. In section III, the coefficients in the effective Lagrangian are determined. 
In section IV we apply the effective three-dimensional field theory to compute the free 
energy in scalar electrodynamics to order A^, Ae^ and e^. In section V we summarize and 
conclude. In two appendices, the notation and conventions are given. We also list the 
sum-integrals used in the full theory as well as the three-dimensional integrals needed in 
the effective theory. 
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2 Effective Field Theory and Dimensional Reduction 



In the imaginary time formalism (ITF) of quantum field theory, there exists a path integral 



representation of the partition function |15|: 



Z = J V(l)exp[- J%T J d^xCs]. (1) 

Here, is a generic field, and £e is the Euclidean Lagrangian which is obtained from the 
usual Lagrangian after Wick rotation, t —it. Moreover, the boundary conditions in 
imaginary time are that bosonic fields are periodic in the time direction with period (3 
and that fermionic fields are antiperiodic with the same period. If $ denotes a bosonic 
field and \1/ denotes a fermionic field we have 

<I>(x, 0) = <l>(x, r), ^(x, 0) = -^(x, r). (2) 

The (anti)periodicity implies that we can decompose the fields into Fourier components 
characterized by their Matsubara frequencies: 

$(x,r) = /3-^5:0„(x)e^'^"^ (3) 

n 

*(x,r) = /?-^^V^„(x)e*'^"^ . (4) 

n 

Here 

ojn = 2mTT, bosons (5) 

ujn = {2n + 1)ttT, fermions, (6) 

where n is an integer. In four dimensions, bosonic fields are assigned a naive scaling 
dimension one, while fermionics field have a naive scaling dimension of 3/2. This im- 
plies that the bosonic fields in the three-dimensional effective theory have a naive scaling 
dimension of 1/2. 

In the ITF, we can associate a free propagator A„(k) with the nth Fourier mode: 



= J^- (7) 



Here, /c = |k|. Hence, a quantum field theory at finite temperature may be viewed as 
an infinite tower of fields in three dimensions, where the Matsubara frequencies act as 
tree-level masses. The n = bosonic mode is called a light or static mode, while the n ^ 
modes as well as the fermionic modes are termed heavy or nonstatic. 

At tree level the light modes are actually massless. If $ denotes the scalar field in 
5f^$^-theory, the scalar field acquires a thermal mass of order gT. The thermal mass 
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reflects that static scalar fields are screened in the plasma, in complete analogy to the 
screening of static electric fields in QED. For the heavy modes this represents a pertur- 
bative correction which is down by a power of the coupling, and these modes are still 
characterized by a mass of order T. However, the light mode is no longer massless, but 
its mass is of order gT. Hence, we conclude that we have two widely separated mass 
scales at high temperature, which are T and gT. The Appelquist-Carrazone decoupling 
theorem [1^| then suggests that the heavy modes decouple on the scale gT. Moreover, 



modern developments in renormalization theory guarantee that static correlators of the 
full theory can be reproduced to any desired accuracy by an effective three-dimensional 
field theory of static mode at long distances R ^ 1/T hj tuning the parameters of the 
effective Lagrangian as functions of temperature and the parameters in the underlying 
theory. The effective Lagrangian is nonrenormalizable and contains infiniely many terms. 
The coupling constants of CqQ encode the physics on the scale T. This process of go- 
ing from a full four dimensional theory to an effective three-dimensional Lagrangian is 
called dimensional reduction. This is the key observation and the starting point for the 
construction of effective field theories at finite temperature. 

Let us now move on to to massless scalar electrodynamics. The Euclidean Lagrangian 
of SQED is 

>^:SQED = \f^uF^, + (P^$)t(I)^$) + ^($t$)2 + £gf + £gh. (8) 

The covariant derivative is = df^+icA^^. In the present work all calculations are carried 
out in the Feynman gauge, but we emphasize that the parameters in the effective theory 
are gauge independent. The gauge fixing term is 

Cgf=\{d,A,)\ (9) 

If f] denotes the ghost field, the ghost term is then 

J^gh=idMd,v), (10) 

and is thus decoupled from the rest of the Lagrangian. 

In 5'^$'^ theory we have seen that we have the two momentum scales T and gT. In 
scalar electrodynamics we have the scale T as well as the scale of electric screening and 
the scale of scalar screening. In the rest of the paper we shall assume that ~ A <^ 1 so 
that counting of powers of the coupling constants is simplified. 

We call the effective three-dimensional field theory electrostatic scalar electrodynamics 
(ESQED), in analogy with the definitions introduced by Braaten and Nieto in the case 
of QCD pT|. The first task is to identify the appropriate fields and the symmetries in 
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ESQED. The fields can be identified (up to normalizations) with the zero-frequency modes 
of the original fields. £esqed contains a real massive scalar field, which we denote by p. 
This field is identified with the n = mode of the timelike component of the gauge field in 
full SQED. There is a complex massive scalar field, 0, which is identifield with $. Finally, 
there is a massless gauge field, which is denoted by A^'^ and is identified with the spatial 
component of the gauge field in the underlying theory. Since the fields of ESQED are 
identified with the static modes of SQED, we have to a first approximation the following 
relations: 



Now, £esqed niust be a gauge invariant function of the spatial fields Aj, up to the 
usual gauge fixing terms. This symmetry follows from the corresponding symmetry in 



the full theory and the Ward-Takahashi identity in the high temperature limit |16|. Since 
SQED is an Abelian gauge theory, there will be no magnetic mass [^. The fact that the 
fields p and $ are massive reflects the screening of static electric and static scalar fields, 
respectively. There is rotational symmetry for p as well a Z2-symmetry for the fields 
and p. The effective Lagrangian then has the general form 

>Cesqed = \F,,F,j + {V4)\V4) + M\K)cj)^ + ]^{dipf + 

\m\K)p' + /i|(A)0Vp2 + ^(0t0)2 + + + (12) 

Here, A is the ultraviolet cutoff in the effective theory. The A-dependence of the param- 
eters is necessary in order to cancel the dependence of A which arises in perturbative 
calculations in electrostatic scalar electrodynamics. The scale A can be thought of as an 
arbitrary factorization scale separating the scale T and eT. Furthermore, 5C represents all 
local terms that can be constructed out of Ai and p, which respect the symmetries of the 
theory. This includes renormalizable terms, such as gE{A.)p^, as well as nonrenormalizable 
ones like hE{A){FijFij)'^. 

It is not too difficult to find the order in the couplings at which a certain operator in 
CqQ starts to contribute to a physical quantity. Take e.g. As(A)p^; The coefficient in front 
of the quartic coupling at leading order in e^, can be found by considering the one- loop 
contribution to the four-point function for timelike photons at zero external momenta. 
Thus, it goes like e^. Moreover, this operator contributes to the free energy first at the 
two-loop level (the double bubble), and each loop is proportional to m£;(A). Since the 
mass goes like eT, this implies that the leading order contribution is of order e^. Other 
operators are analyzed in a similar manner. 

In (|12|), we did not include the unit operator. The coefficient of the unit operator, 
which we denote /esqed(A), gives the contribution to the free energy from the momentum 
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scale T. So if we are interested in calculating the pressure we must determine it, as we 
determine other coefficients in ESQED. Generally, this coefficient also depends on the 
renormalization scale A. By including /esqed(A) in £esqed, we have two equivalent ways 
of writing the partition function in SQED in terms of its path integral representation. In 
the full theory we have 

Z = J VfjVr]VA^V'^Vi;exp[- dr J d^xC\, (13) 
where t] denotes the ghost field. The result using the effective three-dimensional theory is 
Z = e--^ESQED(A)v y VrjVrjVAiVpexp [ - J rf^x^ESQED]- (14) 
The free energy may then be written as the sum of two terms: 



/esqed(A) - IuZesqed/V , (15) 



The ffist term represents the contribution from the scale T and the second term represents 
the contribution from the scale eT, and may be calculated perturbatively using ESQED. 

Let us close this section by pointing out a difference between Abelian and Non-Abelian 
gauge theories. In QCD, we have three momentum scales. The scale T, which is a typical 
momentum of a particle, the scale gT which is the scale of colour electric screening, and 
g'^T, which is the scale of colour magnetic screening. So, it is convenient to construct 
a sequence of two effective field theories in order to unravel the contribution to the free 
energy from the three momentum scales T, gT and g'^T, as demonstrated by Braaten and 



Nieto ||21|. Since the gauge field in ESQED is massless, one could construct a second 
effective field theory only involving Af'^ by integrating out the massive fields p and $. 
However, below we shall argue that this is in fact not necessary. 

We call this effective field theory magnetostatic scalar electrodynamics, (MSQED), 
and it is valid on distance scales R ^ 1/eT. The Lagrangian consists of all terms made 
up by A^'^ which respect gauge invariance and the other symmetries of MSQED. At lead- 
ing order in the couplings of full SQED, MSQED is simply free Maxwell theory in three 
dimensions. However, the only ultraviolet divergences one would encounter doing pertur- 
bative calculations in MSQED are powerlike. The coefficients in front of these depend 
upon the regulator, and are thus artifacts of the regulatorBecause of the unphysical char- 
acter of these divergences, it is convenient to use a regulator where the these divergences 
are subtracted. Dimensional regularization provides such a procedure, since power diver- 
gences by definition are set to zero. This implies that there will be no contribution to the 
free energy from perturbative calculations in MSQED. This fact refiects that there is no 
magnetic screening in Abelian gauge theories and no contribution from the scale e^T to 
the free energy. Hence, there is no need to construct MSQED. 
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3 The Coefficients in the Effective Lagrangian 



In this section we shall determine the parameters of £esqed to the order in A and 
which is required for obtaining the free energy to order A^, Ae^ and e^. We know from the 
previous discussion that static correlators in SQED at long distances R ^ 1/T can be 
reproduced in ESQED by tuning the parameters as functions of A, e^, T and the renormal- 
ization scale A. We determine the coefficients in £esqed by the use of strict perturbation 
theory. This way of determining the parameters of the effective three-dimensional theo- 
ries were first carried out by Braaten and Nieto in Ref. [^|. Strict perturbation theory is 



simply ordinary perturbation theory and is therefore plagued with infrared divergences, 
which become more and more severe as we go to higher orders in the loop expansion. 
Accordingly, we split the Lagrangian of SQED into a free piece and an interacting piece: 

(/:sqed)o = ^F^.F^, + (9^$)t(9^$) + /:gf + £gh, 
(>CsQED)int = e2$t$Aj-zeA^($t9^$-<|.9^$t) + ^($t$)2. (16) 

Although strict perturbation theory breaks down on the scale eT, we can nevertheless use 
it as a tool for determining the parameters of ESQED. The idea is that the coefficient of the 
effective Lagrangian are sensitive to the scale T, but insensitive to the scale eT. We only 
have to be sure that we make the same incorrect assumptions using the effective theory in 
matching calculations. The infrared divergences which occur in perturbative calculations 
in the effective theory will then be exactly the same as those encountered in the full 
theory, and so they cancel against each other in the matching procedure. In ESQED, 
these incorrect assumptions amounts to treating the mass parameters as interactions. 
The partition of the Lagrangian of ESQED is 

(>Cesqed)o = ^FijFij + {di(j))\di(p) + ]^{dipf + ^gf + ^g^, 
(^ESQED)int = ^m|(A)/ + M2(A)0V + e|(A)0t0Af + 4(A)0Vp2 



-zes(A)Af (0^9,0 - <pd4^) + ^^(0V)' + 5C. (17) 

Since the strict pertubative expansion given by ( [16|) and ( |T7| ) is naive perturbation theory, 
it should be clear that the coupling constants of ESQED can be written as power series 
in A and e^. 

In full SQED, wiggly, dashes and dotted lines denote the propagators of photons, 
charged scalars and ghosts, respectively. In ESQED, the same conventions apply. More- 
over, solid lines denote the propagators of the real scalar field p. 
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3.1 Coupling Constants 



In the present work, one needs the couphng constants A3 (A), e|;(A) and h'^{A) to leading 
order in A and e^. At this order one can simply read off these coefficients from the 
Lagrangian in full SQED. We have previously found that p = -^^o , = "^"^ 
Af'^ = -^Ai at leading order. Substituting these expressions into (H) and comparing it 
with (|1^), we conclude that 

A3(A) = AT. (18) 

Similarly, one finds 

e|(A) = e^T, 4(A) = e'T. (19) 
There is no dependence on the renormalization scale A. 

The above parameters have been obtained at next-to-leading order in A and in e.g. 
Ref. p6l by considering the appropriate static correlators. 



3.2 Mass parameters 

We shall now determine the mass parameters m|;(A) and M^(A) appearing in the effec- 
tive Lagrangian. These parameters are needed to leading order in A and e^. The mass 
parameter m|;(A) has previously been calculated by the present author in the two- loop 
approximation using strict perturbation theory [^. The mass parameter M^(A) has also 
been calculated to two-loop order in Ref. using the effective potential. Nevertheless, 



we include the calculations here for completeness. 

The mass parameters have the following interpretation; they represent the contri- 
bution from the scale T to the electric and scalar screening mass, respectively. These 
screening masses give information of the screening of static electric and scalar fields in 



the plasma Now, the scalar screenning mass cannot be calculated beyond leading 
order due to a mass shell singularity, as discussed by Blaizot and lancu in Ref. The 
problem is the same as in QCD, namely a scalar field interacting with a massless gauge 
field in three dimensions. In QCD, this singularity may be screened by a magnetic mass 
of nonperturbative origin. In SQED the magnetic mass is absent since it is an Abelian 
theory |2^, and so the problem cannot be solved this way. We shall not discuss this any 
further, but refer to Ref. |^ where a nonperturbative definition of the scalar screening 
mass is discussed in detail. Nevertheless, these infrared problems do not prevent us from 
determining the mass parameter M^(A) using the strict perturbation expansion. 

There are several ways of determining the mass parameters, and the simplest way to 
do so is to match the screening masses in the two theories. The screening mass is defined 
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as the pole position of the propagator at spacehke momentum pO|. Let us denote the 



static self-energy function of the zeroth component of the gauge field by n(k). Similarly, 
we denote the static scalar self-energy function by S(k). The corresponding screening 
masses are and respectively. In the full theory, one then finds that the serening 
masses are solutions of the following set of equations: 

+ n(k) =0, = -ml, (20) 

fc2 + S(k) = 0, k^ = -Ml (21) 

The self-energy functions can be expanded in powers of the external momentum k and 
also in the number of loops in the loop expansion. The nth order contribution to the self- 
enregy functions are then denoted by 11*^"') (k) and E*^")(k), respectively. At leading order 
in the couplings, we simply substitute the one-loop approximations to the self-energy 



functions into (pOD and (^). The solutions are 

m^^n«(0), M2^S«(0). (22) 

The symbol ^ is a reminder that the expressions are obtained using strict perturbation 
theory. Since strict perturbation theory treat infrared effects incorrectly, the mass param- 
eters generally do not coincide with the physical screening masses . At leading order, 
however, they normally do coincide. 

In the full theory the Feynman graphs in Figs. |l| and ^ display the contributions at 
the one-loop level to the the self-energy function of p and $, respectively. We find 

n'"(0) « 2e'f^±-4e^f^^,, (23) 

E<"(0) « id-De'f^-L + fgj,. (24) 

Here, d = 4— 2e. In ESQED, the self-energy function of the real scalar field p is denoted by 
ngff(k. A), while the corresponding self-energy function of (j) is Sgff(k, A). In the effective 
theory the screening masses are by the following set of equations, where the superscript 
"(1)" indicates that the self-energies are taken in the one-loop approximation: 

k^ + m\A) + U^^{0,A) = 0, k^ = -ml (25) 
e + M\K) + S^'^(0, A) = 0, k^ = -Ml (26) 

In the effective theory, all diagrams vansih identically; Since all fields are massless in strict 
perturbation theory, there is no scale in the loop integrals when the external momentum 
k is set to zero. The integrals are set to zero in dimensional regularization. Hence, 
Hgg(0,A) = and Sg^(0,A) = The matching conditions then imply that m\{K) = 
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n(i)(0) and M2(A) = S(^)(0). Using (^)-(^]5D of Appendix A, we obtain the mass 
parameters at leading order in A and e^: 



m|(A) 



+ 



AT2 



4 18 

To leading order in the couplings, then, the mass parameters are independent of A. 



(27) 
(28) 



3.3 Coefficient of the Unit Operator 

In this subsection we shall determine the coefficient of the unit operator to next-to-leading 
order in A and e^. We shall consider the one, two and three- loop contributions separately. 
The matching condition we use to determine /esqed(A) follows from the two path integral 
representations of the partition function 

\nZ = -/esqed(A)1^ + In^ESQED. (29) 

Let us first focus on the SQED. The one-loop graphs in the underlying theory are depicted 
in Fig. ^ and the corresponding contribution reads 

This is simply the free energy of a noninteracting gas of photons and charged scalars. 

The two-loop diagrams are shown Fig. ^ After some purely algebraic manipulations, 
they factorize into products of simpler one-loop sum-integrals. Using Appendix A, we 
obtain 

^^^^ 1 3 1 r^'^'r^-L^e^ 

^fpg + - 2^^^'' tpQ = ( 12 ) is + • 

Here, Z\ and Z^.^ are the renormalization constants for the quartic coupling and the gauge 
coupling, respectively. 

Let us now turn to the three-loop diagrams. These are displayed in Fig. |^. Here, 
the shaded blob means insertion of the one-loop photon polarization tensor H^,^(/co,k), 
while the black blob implies insertion of the scalar self-energy function E(/co,k), also in 
the one-loop approximation. 

The first four diagrams can expressed entirely in terms of the bosonic basketball, 
which is defined in Appendix A. After some purely algebraic manipulations involving 
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several changes of variables, one finds: 



rA 



The fifth diagram gives a contribution 



PQK P^Q^K^{P + Q + K)^' 



1 

4 Tpqk P^Q^K^{P + Q + Ky 



(32) 



Tpqk P^Q^K^{P - Ky{P - Qf 
Tp, 



fpQK P^Q'^K^ ' 



(33) 



The last graph reads 
1^ 1 



2rp P4 



S(po, P) 



-2A2 



L 9 



+ 



2(rf-l)Ae2 (rf-l)V- 



-2e 



+ 
1 



PQK P^Q^K^ 



^4 . 

Tpqk P^Q^K^P + Q + Kf 



(34) 



Note that the fifth as well as the sixth diagram are infrared divergent. It stems from the 
facts that both the zeroth component of the static photon polarization tensor, noo(0,k), 
and the static scalar self-energy function, S(0, k), are nonzero in the infrared limit k — > 0. 
These diagrams are the first in an infinite series of infrared graphs which are called ring 
diagrams or plasmon diagrams. They are obtained from the two diagrams considered 
above, by one or more insertions of T{^y{kQ, k) and S(/i;o, k), respectively. 

Renormalization of the quartic coupling is carried out by substituting Z\ into (|3T|). 
To order A^, Ae^ and e^, the renormalization constant reads 



ZaA = A + 



5A^ - 18Ae^ + 54e^ 
487r2e 



(35) 



Similarly, we renormalize the gauge coupling by subtituting Z^a into (p]). To order e^, 
Zf.2e^ is given by 

(36) 



487r2e' 

The result in full SQED is now given by the sum of ^ After we have carried out 
coupling constant renormalization, we are still left with terms proportional to 1/e: 



TlnZ 
V 



45 



+ 







5e2n 




[3^ 


^^1 



13 



A2 ,T\2rl(} 



167r2 V12/ 



A 



In 

9 AttT 



167r2 Vl2y 

—i-f 

167r2 V12/ 
436C'(-i; 



4 31 2 C^(-3) 4 CX-1) 
^9^^ ^ 45 3 C(-3) 3C(-1) 
A 



- + 201n^ + 4,. + |- 

18 365 , A 

— + — In -— + 137i, + 
e 3 AttT 



16 



C(-i) 

236 110C'(-3) 
~3 3~C(-3) 



(37) 



Again, the sign 



3 C(-l, 

indicates that the resuh is obtained using strict perturbation theory. 



Let us next consider ESQED. The contributions to Iu^esqed comes from ordinary 
one. two and three-loop diagrams, as well as one-loop graphs with one or two mass inser- 
tions. The loop-integrals all vanish in the strict perturbation expansion, since there is 
no mass scale. However, note that some of the three-loop diagrams as well the one-loop 
graphs with two mass insertions and the two-loop graphs with a single mass insertion are 
linearly divergent in the infrared. These divergences are identical to those encountered in 
full SQED (the fifth and sixth diagrams in Fig. |^). 

The vanishing of Iu^esqed implies that the free energy in the strict perturbation 
expansion is given by 

TlnZ 



V 



T 



/esqed(A) + (5/esqed(A) 



(38) 



Here, 5/esqed(A) is the counterterm of the unit operator. According to 5/esqed(A) 
can be determined by calculating the logarithmic ultraviolet divergences in the effective 
theory. Generally, 5/esqed(A) is given by a polynomial in the parameters of ESQED. At 
leading order it turns out that it is given by 



(^/esqed(A) 



e|M2 1 
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2(47r)2 e' 

which follows from a two-loop calculation in the next section. Since the mass M is 
multiplied by 1/e, it is necessary to expand it to first order in e when expressing ^/esqed(A) 
in terms of e^, A and T. From (E3D one finds 
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Putting our results together, one finally obtains 
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We have here used the renormalization group equations 

dX 

to change the scale from A to /i. 



247r2' 

5A^ - ISAe^ + 54e^ 
24^^2 ■ 



(42) 

(43) 
(44) 



The remaining logarithm of A shows that /esqed(A) 
depends explicitly on A. This logarithmic dependence on A is necessary to cancel the A- 
dependence which arises in calculations in ESQED, as we shall see in section IV. This is 
contrast with spinor QED , where the unit operator is independet of the factorization 
scale to order e"^. 



3.4 Renormalization Group Equations in ESQED 

Above, we have seen that some of the parameters in ESQED do not depend explicitly 
on the cutoff A to the order which we calculated them, while others do. More generally, 
the coefficients in £esqed satisfy a set of renormalization group equations. The solutions 
to these equations can be used to sum up leading logarithms of the couplings constants 
from higher orders in the perturbation expansion, as discussed in detail in Ref. . 
Although we shall not need these evolution equations in the present work, we include a 
short discussion for completeness. They will be relevant when calculations are pushed to 
higher orders. 

The renormalization group equations follow from the requirement that physical quan- 
tities such as the electric screening mass and the free energy be independent of this 
arbitrary factorization scale. If C„(A) denotes an arbitrary operator, one can write 

A^^ = PniCiA)). (45) 
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The beta functions /3„(C(A)) can written as power series in the couphng constants of 
the effective theory. Using dimensional arguments we can infer the general structure. 
Consider first the coefficient of the unit operator /esqed(A) and let us restrict ourselves 
to the superrenormalizable terms which are explicitly displayed in (|l2l). The unit operator 
has dimension three. The beta function then contains terms in the form e\m\, h\m\, 
\^m\, e|;M^, h^M"^ and AsM^, as well a cubic polynomial in the three coupling constants 
e%, h\ and A3. To order Ae^, Ae^ and only the first kind of terms can appear. The 
coefficients in front of these terms are determined by calculating the logarithmic ultraviolet 
divergences in ESQED. This is done in the next section, and the result is 

^dks^^-ijM^^O(e% (46) 

Using the fact that e\ = e^T at leading order and the leading order result for M^, we can 
write the equation for /esqed(A) in terms of A and e^: 



feQED(A) _ .e^ Ae^ 
dk (47r)2^2 9 



^u^j^^^^^ = -!-[;- + ^] + 0{e% (47) 



This form of the renormalization group equation could of course have been obtained 
simply by differentiating (|42D with respect to A. 



Consider next the mass parameter m\{A) and M^(A). The beta functions for these 
parameters must be quadratic polynomials in the A3, e\ , h\ and A^; (the latter is the 
coefficient in front of p^). We have already noted their independence of A at leading order 
in the couplings. In Ref. |2^, I have calculated m\[K) to next-to-leading in A and e^. 
The mass parameter turns out to be independent of the factorization scale also to this 
order in the couplings: 

A^4^ = 0{e% (48) 

This is in contrast with the corresponding result for the scalar mass parameter M^(A). 
In Ref. ||2^ , this parameter has also been calculated at next-to- leading order in A and , 



and the beta function reads 



dM\K)_ _ _J_ ^^^4 4A3e| , 2\l 
dA ~ jAn) 



A^^^ = 7^[6e^ - ^ + ^] + Oie% (49) 



Finally, we make the remark that the f3 functions for the couplings A3 (A), e|;(A), h'^{A) 



and Aa(A) vanish to order A^, Ae^ and [^. This property actually holds to all orders 
in superrenormalizable interactions. 

In summary, there are two parameters in ESQED whose evolution are relevant up to 
order e^. These are /esqed(A) and M^(A). The solutions to the renormahzation group 
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equations (|46|) and ( ^9]) are, respectively 

/esqed(A) = /ESQED(A')-^|^ln^, (50) 
M\A) = M^(A0+(^[6e|-^ + ^]ln|,. (51) 

4 Calculations in ESQED 

Having calculated the coefficients in ESQED to the order required, we are now ready to 
calculate the free energy in SQED to order A^, Ae^ and e^. The free energy is given by 
the two terms in ([I5|), where /esqed(A)T is given by (^21). The remaining term IuZesqed 
is given by perturbative calculations using £esqed- In order to obtain this term, it is 
necessary that we take screening effects properly into account. This is done by including 
the mass parameters m|;(A) and M^(A) into the free part of the Lagrangian. We then 
include the effects of the mass parameters to all orders, while all other operators are 
treated as perturbations. Thus, we split ^esqed- according to: 

(/:esqed)o = ^F,,F,^ + (9,0t)(5,0) + M2(A)0V+l(a,p)2 

+M^(0t. 



The one and two-loop contributions to the pressure are depicted in Fig. ||. By adding 
the counterterm 5/esqed(A), we obtain the contribution to the free energy from the 
momentum scale eT. 



(>CESQED)int = e|(A)0t0Af + ^l(A)0Vp' - *e^(A)^f (</'^9,</> - 09, 

+ 6C. (52) 



T In 2^ESQED 
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It f ln(p2 + rnl)+T f ln{p^ + M^) + - 2)T / \np^ 
z Jp Jp Z Jp 



2 Jpg {p^ + M^){q^ + M^){p - qy ^ Jpg {p^ + M^)q^ 

9rr. f 1 r 1 

+e|T / ^^^^^..^^^^ + 4t 



pg (p2 + M2)(g2 + m|) 3 Jpg (p2 + M2)(g2 + M2) 



+T5/esqed(A). (53) 

Here, c? = 3 — 2e. The two-loop contributions may be reduced to products of one-loop 
integrals and a two-loop integral, which are tabulated in Appendix B. This yields 

TIuZesQED _ _^_M!l^^lMl^^lMl\l^O^A^ Al^ 2^^ 

V ~ 127r 67r + 327r2 + 327r2 ^e^^ 2mJ + ''^ 167r2 
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+ 



487r2 



+ 5/esqed(A). 



(54) 



The logarithmic ultraviolet divergence stems from the theta-diagram. This pole in e is 
then canceled by the counterterm of the unit operator, 5/esqed(A), which is given by (|39|). 
Using the expressions for the the parameters in ESQED in terms of the couplings in full 
SQED as well as the factorization scale A and temperature, we get the contribution to 
the free energy from the scale eT: 
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(55) 



Adding (^) and (pSf), we finally obtain the free energy of massless scalar electrodynamics 
at high temperature through order A^, Ae^ and e^: 
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(56) 



This represents the main result of the present paper, and at this point some comments 
are in order. Firstly, we note that the two-loop contribution in the gauge sector is exactly 
the same as the one appearing in spinor QED |jl3[. The terms which go like and A^^^ 
are non-analytic in the coupling constants, and therefore demonstrate that one receives 
contributions from all order in naive perturbation theory. The above terms represent the 
leading order contributions from the infinite string of infrared divergent ring diagrams. 
We also note that the A-dependence cancels to next-to- leading order in A and e^, leaving a 
result which is renormalization group invariant to order A^, Ae^ and e^. This can be easily 
checked by using the renormalization group equations, ( ^5] ) and (^^, for the running 
coupling constants. 

The effective field theory approach also explains the occurrence of logarithms of T/M 
in the expression for the free energy. We have seen that /esqed(A) depends explicitly 
upon the factorization scale A and such terms appear as logarithms of A/T. In ESQED, 
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logarithms in the form A/M occur in perturbative calculations. In order to cancel the 
dependence of A in the final answer, these logarithms must match the logarithms from 
/esqed(A). We are then left with logarithms of M/T. Hence these logarithms are con- 
nected to the renormalization of /esqed(A). 

A similar term g'^hi{mE/T) has also been also found by Braaten and Nieto in the 
case of QCD No term of order g^\\i{m/T) arises in gf^^^-theory, since /(A) does not 



run at next-to-leading order in g |20|. This remark also applies to QED p 



5 Concluding Remarks 

In the present work we have calculated the free energy in SQED to order A^, Ae^ and 
by effective field theory methods. The major advantage of the effective field theory 
approach is that it allows one to work with a single scale at a time. This makes per- 
turbative calculations simpler, since the scales T and eT are not intertwined in already 
complicated sum-integrals in the full theory. However, we would also like to emphasize 
that effective field theory have some limitations. At present it is not clear how to apply it 
in the computations of dynamical quantities. Here, the resummation program of Braaten 
and Pisarski or extensions thereof are still essential ingredients of high temperature field 
theory. 

We would also like to outline the calculations necessary to push the calculations of 
the free energy to one more order in the coupling constants. The coefficient of the unit 
operator are written as a power series in A and e^, and so we already know this parameter 
to sufficient precision. The free energy at one-loop in the effective theory goes like M^(A)T 
and m\{A)T. This implies that we need to know the mass parameters in ESQED at next- 
to-leading order. The former has been calculated by Farakos et al. in Ref. p6[ while the 



present author has computed the latter in Ref. p2|. Furthermore, there are no new 
operators which contribute at this order, so all that remains is calculating the three-loop 
diagrams in ESQED. As we have already noted both the coefficient of the unit operator 
and the scalar mass parameter depend explicitly on the scale A. We realize that, upon 
expanding M^(A) in powers of coupling constants, we get terms proportional to In A/T. 
This allows us to predict terms in the form In M/T or In m^j/T in the pressure of 
scalar electrodynamics. This also explains the absence of such terms in QCD I^T] as well 



as QED since the corresponding mass parameter m%{K) is independent of A to fourth 
order in the gauge coupling. 

Finally, we would like to briefly comment upon other applications of the present 
approach. This method is perhaps the most transparent way of doing effective field theory 
at finite temperature, and it would certainly be worthwhile extending it to e.g. the study 
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of phase transtions. This aproach would then represent an alternative to the methods 
which explicitly construct the effective field theory of the light mode by integrating over 
the heavy modes . Work along these lines is in progress . 
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A Sum- integrals in the Full Theory 



In this appendix we give the necessary details for the sum-integrals used in present work. 
Throughout the work we use the imaginary time formalism, where the four-momentum is 
P = {po, p) with P"^ = pI + p"^. The Euclidean energy takes on discrete values, po = 2n'7iT 
for bosons. Dimensional regularization is used to regularize both infrared and ultraviolet 
divergences by working in (i = 4 — 2e dimensions, and we apply the MS renormalization 
scheme. We use the following shorthand notation for the sum-integrals that appear below: 



P0=2iTnT 



(2vr; 



3-2e 



fiP)- 



(A.l) 



Then fi coincides with the renormalization scale in the MS renormalization scheme. 



Arnold and Zhai |jT2| have developed a new machinery to deal with complicated sum- 
integrals. The methods are completely analytic and represent significant progress in 
perturbative calculations at finite temperature. They have calculated and listed all the 
sum- integrals needed in the present work, except for the one in ( |A.8| ). This sum- integral 
has been evaluated by Braaten and Nieto using the methods of Ref. . We repro- 
duce them below for the convenience of the reader. 



The specific one-loop sum-integrals needed are 
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All of the two-loop sum-integrals encountered in this work factorize into products of the 
one-loop sum-integrals listed above. The only two-loop sum-integral needed is: 



1 



PQ P^Q\P + QY 



0. 



(A.6) 



Some of the three-loop sum-integrals factorize into products of three one-loop sum- 
integrals, while others factorize into a product of two one-loop sum-integrals and the 
two-loop sum- integral in ([A.6|) . The specific three- loop sum-integrals needed in the present 
work are 
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The first of these three-loop sum- integrals is the bosonic basketball 1T2|, while the second 
is denoted by 2, -2 in Ref pT| . 



B Integrals in the Effective Theory 

In the effective three-dimensional theory we use dimensional regularization in 3 — 2e 
dimensions to regularize infrared and ultraviolet divergences. In analogy with Appendix 
A, we define 

Again coincides with the renormalization scale in the modified minimal subtraction 
renormalization scheme. 

All the integrals needed in the present work have been evaluated and listed by Braaten 



and Nieto in Ref. 20 . We list them below for the convenience of the reader. 
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In the effective theory we need the following one-loop integrals 



/ ln(p2 
Jp 



2m 3> 



= -^[l+(21n-^ + 2)e + 0(e2)l, (B.3) 



/p + m^ 

f T^T^—^ = + (21n7^)e + 0(e')l. (B.4) 

Jp{p^ + m?Y SvrmL V 2m^ ^ ^ ^ 

The two-loop integrals in the effective theory also factorize into products of one-loop 
integrals plus a single two- loop integral. More, specifically one needs 

1 1 r 1 1 , /i X 

+ - + ln-^ + e 
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Figure 1: One- loop scalar self-energy diagrams in SQED. 




Figure 2: One-loop self-energy diagrams in the full theory. 




Figure 3: One-loop vacuum diagrams in SQED. 
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Figure 5: Three- loop diagrams contributing to the free energy in the underlying theory. 
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